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Abstract 

The theory of strong interactions, quantum chromodynamics (QCD), is quite successful in the 
prediction and description of main features of multiparticle production processes at high ener- 
. . . gies. The general perturbative QCD approach to these processes (mainly to e^e~ -annihilation) is 

. briefly formulated and its problems are discussed. It is shown that the analytical calculations at 

\ the parton level with the low-momentum cut-off reproduce experimental data on the hadronic final 

CN ' state in multiparticle production processes at high energies surprisingly accurately even though the 

perturbative expansion parameter is not very small. Moreover, it is important that the pertur- 
bative QCD has been able not only to describe the existing data but also to predict many bright 
qualitatively new phenomena. 
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1 Introduction 

Multiparticle production is the main process of very high energy particle interactions. Studying it, 
one hopes to get knowledge on the validity of our general ideas about the structure of the matter at 
smallest distances, on new states of matter which could be created at these extreme conditions, on 
the asymptotic properties of strong interactions, on the confinement of quarks inside hadrons etc. One 
should also understand these processes to separate the signals for new physics from the conventional 
background, in particular, from the features explainable in the framework of quantum chromodynamics. 
This is necessary for the successful planning of new experiments as well. 

With studies of cosmic rays and the subsequent steady increase of energies of particle accelerators, the 
interest to multiparticle production processes increased as well. The theoretical interpretation of these 
processes passed through several stages, in detail described in many monographs and review articles. It 
evolved from the statistical and hydrodynamical approaches to the peripheral, multiperipheral, Regge, 
parton models and QCD. Let us stress that all of them are somehow used now (quite often in various 
combinations for the description of different stages of the process) in the computerized Monte Carlo 



The ideas of quantum chromodynamics, where quarks and gluons play a role of partons, nowadays 
prevail even in the phenomenological models. Neither these phenomenological models nor the widely 
used analytical perturbative QCD approach can evade an important problem of the transition from 
partons (quarks, gluons) to the observed particles. This stage is treated phenomenologically within 
some hypotheses in any of them. It introduces additional parameters which, in principle, can give us a 
hint to the confinement property, but they sometimes are hard to control and directly extract. 

The whole process of multiparticle production is considered as a jet-like cascade of consecutive emis- 
sions of partons each of which produces the hadrons observed. Jets from primary quarks were discovered 
in e"'"e~-collisions in 1975 with the angular distribution expected for spin 1/2 quarks. Gluons emitted 
by quarks at large transverse momenta can be described by perturbative QCD due to the asymptotic 
freedom property according to which the coupling strength in QCD decreases with the transferred mo- 
mentum increasing. Such processes are used to determine the value of the coupling strength. However, 
one can try to proceed to lower transverse momenta when many jets (and consequently many hadrons) 
are created. These processes are of main concern for this survey. 

Main attention will be paid to the process of e"'"e~-annihilation at high energies, where two point- 
like particles collide and the initially created state is determined by the time-like quark- ant iquark pair 
whose evolution gives rise to the production of additional jets. Perturbative QCD has been widely 
used for the treatment of these processes (as well as for hard jets created in the final states of ep, pp 
and similar interactions). Its predictions and comparison with experimental data are described in the 
monographs |l], 0, 0, ^, |^ and recent review papers 0, |^, ||, |, |10[- The application of the QCD ideas 
to comparatively soft processes in hadronic and nuclear collisions requires for some modification of the 
approach with account for the internal structure of colliding objects (structure functions) as this is 



described, e.g., in the surveys |TT|, |12|, |13|, |14|, 0, |l| 



To be compared with experimental data, the results of quantum chromodynamics should be pre- 
sented either by analytical formulas or numerically as computer calculations according to the Monte 
Carlo models. 

Usually, the Monte Carlo models deal with matrix elements (actually, with the probabilities) of 
a process at the parton level plus the hadronization stage which uses either the string model |T] 



or the fragmentation of hadronic clusters [|T^, [T^. They properly account for the energy- momentum 



conservation laws because all characteristics of the exclusive process are calculated by computing. One 
gets all possible characteristics at a given energy but fails to learn their asymptotic behaviour. Let 
us note that many Monte Carlo models have been proposed which differ by the relative role of the 
parton and hadronization stages and, correspondingly, by the set of adjusted parameters, which can be 
interrelated and sometimes are hard to control. Some freedom in their choice also defines the difference 
in their predictions. 

On the contrary, the analytical approach in QCD pretends to start with asymptotic values, where the 
energy-momentum conservation laws can be neglected in the first approximation, and proceeds to lower 
energies accounting for conservation laws, higher order perturbative and simplified non-perturbative 
effects. In the analytical calculations, the perturbative evolution of the jet with the virtuality of the 
initial parton degrading due to emission of secondary partons is terminated at some low scale Qq of 
the order of some hundreds MeV for transverse momenta or virtualities of partons. Some observed 
variables, e.g., such as the thrust or energy flows, are insensitive to this "infrared" cut-off, i.e. they do 
not change if collinear or soft partons are additionally emitted. 

For others, like the inclusive distributions, the local parton-hadron duality (LPHD) is assumed EC 



which declares that the calculated distributions at the parton level describe the hadron observables 
up to some constant factor. This concept originates [^, ^ ^ from the preconfinement property, 
i.e. from the local compensation of colour charges, and, consequently, from the tendency of quarks 
and gluons to form the colourless clusters. Surely, the similarity of the distributions of partons and 
hadrons substantially owes to the fact that physics of the colour confinement is governed by rather soft 
processes with small momentum transfers. This works surprisingly well when applied for comparison 



of the theoretical parton distributions with the hadron characteristics observed in experiment. No 
phenomenological hadronization model is used in this case. At the same time, some specific effects, e.g., 
the Bose-Einstein correlations for the identical mesons, are out of the scope of this treatment. This is 
only admissible if these effects are reasonably small. 

Even more amazing look two other features of the perturbative QCD approach: the probabilistic 
description and its applicability to the comparatively soft processes with low transverse momenta^. In 
the physical gauge, the leading (at high energy) terms appear as the squared moduli of the amplitudes, 
i.e. as the probabilities. At high energies, the quantum interference of different amplitudes for the parton 
production results in the angular (or, more precisely, transverse momentum) ordering of successive 
emissions of gluons which, in its turn, does not spoil the probabilistic equations for these processes. 

The solutions of these equations obtained via the modified perturbative expansion sometimes seem 
to be applicable even for rather soft processes where the perturbative expansion parameter is not 
sufficiently small and, moreover, it is multiplied by some large factors increasing with energy. Thus, 
such an approach can be justified only because some subseries of the purely perturbative expansion 
ordered according to their high energy behaviour are summed first (from here originates the term 
"modified perturbation theory") and then the asymptotic series is cut off at the proper order. 

Surely, the probabilistic treatment is just the simplest approach to the whole description of the 
process valid where the colour coherence is not important, the non-linearity of the process, induced, 
e.g., by the high density of gluons, and, consequently, the unitarization are not taken into account etc. 
If one has to take these effects into account, then the direct consideration of the traditional Feynman 
diagrams is preferred because the influence of these speciflc effects in terms of the generating functional 
technique is hard to treat. Unfortunately, in this case it becomes more difficult to proceed to the 
higher order approximations as well. In this framework of the combined study of the equations for the 
generating functionals and of the interference effects implied by Feynman graphs the perturbative QCD 
has demonstrated its very high predictive power. 

In this paper, we briefly describe the main obtained physical results about the processes of the 
production of hadronic jets, mostly in e"'"e~-annihilation at high energies, without the detailed discussion 
of any speciflc particular technicalities both of the theory and of the experiment. The main goal is to 
provide the reader with a short guide in this vast region without omission, if this is possible, of main 
important results and no hiding of newly appearing problems. I am deeply sorry if something essential 
has been omitted and ask for an excuse from those authors whose papers did not flt the limited space 
of this review paper. 

For example, what concerns the theoretical approaches, after their brief description we demonstrate 
only the main equations and their solutions with more detailed discussion of the obtained results, 
their correspondence to experimental data and physical meaning. This would allow the reader to get 
quickly acquainted with the general situation and then, following the numerous cited original papers, 
more voluminous surveys and monographs, to study it in more detail if some interest to the particular 
problem has appeared. The comparison of theoretical results with experimental data is demonstrated 
in some flgures which constitute only a small part of those available in the papers referred to. 

Thus this review paper is mainly aimed at those who want to learn about the achievements and 
problems in the applications of quantum chromodynamics to the description of multiparticle production 
processes at high energies, especially, those who start their study of high energy particle physics or 
specialists in other flelds. 



2 QCD equations 



The most general calculational approach to the characteristics of multiparticle production processes 
starts from the expression for the generating functional [p4|, p5|, the equations governing its evolution 



^For more details see Refs. [5, 9| 



in the framework of the considered theoretical scheme and the subsequent solution of these equations. 
The generating functional contains the complete information about any multiparticle process and is 
defined as 

G{u,y) = Y^ / (fki...(fknu{ki)...u{kn)Pn{ki,--,kn;y), (1) 

where Pn{ki, fc„; y) is the probability density for the exclusive production of n particles with momenta 
ki, kn at the initial virtuality (energy), proportional to exp(y), and u{k) is an auxiliary function. For 
u{k) =const, one gets the generating function of the multiplicity distribution Pn{y)'- 

n=oo 

Giu,y)= Y: u'^Pniy). (2) 

n=0 

The variational derivatives of G{{u}) over u{k) (or usual u-derivatives for constant u) provide any 
inclusive and exclusive distributions (in particular, the average values and correlators of any rank), 
i.e. the complete information about all characteristics of the process. This information is, however, 
unavailable unless the equations for the evolution of the generating functional are formulated. These 
equations depend, naturally, on the theory which has been used for the description of these processes. 

Two equations have been proposed in quantum chromodynamics to describe the multiparticle pro- 
duction processes in the different kinematical regions. In brief they are named by the first letters of the 
names of their authors. They are the DGLAP equation ^ |28| and the BFKL equation [^, ^ |3T| . 



The DGLAP equation applies to the process of the evolution of the virtualities (or transverse momenta) 
of the system. The BFKL equation deals with this evolution in terms of the longitudinal momenta. 
Correspondingly, the DGLAP equations are used for the description of the development of the partonic 
jets where the initial highly virtual time-like partonQ evolves, e.g., as in e"'"e~-annihilation at high ener- 
gies or hard jets in hadronic processes (see the review papers |P, |^, H, |^). In its turn, the longitudinal 
momentum evolution is typical for the softQ exchange in the t-channel of the interacting particles which 
is of the multiperipheral character with some rescatterings (see the review papers |1TT| , |12|, |13|, |14|, |15], 0). 



Some proposals for the interpolation between these two regions have been promoted |32, 33 . 

Initially, these equations were derived in the framework of the Feynman diagram technique by the 
summation of certain subseries from the whole perturbation series expansion. Their validity has been 
proven in the leading and next-to-leading approximations of the modified perturbation theory (in more 
detail it will be described later). Both equations are applicable to the processes in which the density of 
the gluon field is comparatively low and the multiparton interactions and screening may be neglected. 
They have been successfully applied to interpretation of many experimental results and for prediction 
of new properties of various distributions of secondary particles. 

However, these simplified versions of both equations have many unsolved problems. They are partly 
described below in the case of the DGLAP equation. What concerns the BFKL equation, let us just 
mention only one, however important, problem related to the prediction of the power-like increase 
of the total cross sections with energy which violates the unitarity condition and contradicts to the 
Froissart bound. According to this bound, the total cross sections of hadronic processes are allowed to 
increase with energy s not faster than a oc log^ s. In the leading order of the modified perturbation 
theory, the BFKL equation predicts the power-like increase of the total cross sections with very high 
exponent a oc s°'^, that strongly contradicts not only to theoretical restrictions but also to all modern 
experimental results. Recently, it was shown |3^ that the corrections due to the next order terms 
are quite large, and even though the power dependence of cross sections on energy persists, the exponent 
becomes much lower cr oc s^'^"^ . Such an increase is already much closer to the phenomeno logical fits of 
experimental data at presently available energies which give rise to the exponent values in the range 
from 0.08 to 0.12. One can hope that the higher order corrections will lead to further decrease of the 



•^With large positive 4-momentum squared > 0. 

•^When the longitudinal momentum transfers x are shared by small portions; a; ^ 1. 



exponent while the summation of the subseries of the perturbative expansion would be able to avoid 
the contradiction with the Froissart condition. 

The unified field-theoretical description of both approaches can be obtained if the proper effective 
action of the theory ^ ^ ^ is constructed. In particular, in the papers ^7\, |3^, |39 



apart from the term of the standard Yang-Mills gluonic action, the effective action contains the term 
corresponding to the non-abelian eikonal interaction of the fastly moving sources with the relevant 
component of the gluon field. Such field-theoretical scheme with the effective action admits the gen- 
eralization of these equations to the non-linear case. In the framework of the Wilson renormalization 
group approach one can derive the unified functional QCD equation |37, ^ in the leading logarithmic 
approximation which takes into account the higher orders of the parton density fluctuations. It leads 
to the system of interconnected equations for the partonic correlators of an arbitrary order. However, 
a lot of work is still needed in this direction as to develop the higher order approximation approach as 
to understand the role of the mechanism of the non-linear contributions and the unitarity bounds. 

As previously mentioned, this review is mainly devoted to the e"'"e~-annihilation processes where the 
experimental data are most precise and extensive ones. Therefore the solution of the DGLAP equations 
for the evolution of partonic jets constitutes the main content of the theoretical approach adopted hereQ. 

The general structure of the equation for the generating functional in QCD describing the jet evo- 
lution for single species partons can be written symbolically as 

G" ~ y asK[G ® G - G]dQ. (3) 

It shows that the evolution of the functional G indicated by its derivative G' over the evolution parameter 
(the transverse momentum or the virtuality) is determined by the cascade process of the production of 
two partons by a highly virtual time-like parton (the term G ® G) which provides new partons in the 
phase space volume considered dfl and by the escape of a single parton (G) from a given phase space 
region. 

Therefore this equation contains terms corresponding to inflow and outflow of partons. In fact, it 
can be interpreted as the kinetic equation with the collision integral in the right hand side. The weight 
factors are determined by the coupling strength as and the splitting function K which is defined by 
the interaction Lagrangian. The integral runs over all internal variables, and the symbol ® shows that 
the two created partons share the momentum of their parent. The initial condition for equation (^) is 
defined by the requirement for the jet to be created by a single initial parton, i.e., by 

Pn = Sni; Go = u{k). (4) 

It is clear from this formula that we have to deal with the non-linear integro-differential probabilistic 
equation with shifted arguments in the G ^ G term under the integral sign. 

For quark and gluon jets, one writes down the system of two coupled equations. Their solutions give 
all characteristics of quark and gluon jets and allow for the comparison with experiment to be done. 
Let us write now them down explicitly for the generating functions: 

G'a = f'dxKSixh',[GGiy + \nx)GGiy + Hl-x))-GGiy)] 
Jo 

+ Uf f dxK^{x)^l[GF{y + \nx)Gp{y + ln(l - x)) - Gg(i/)], (5) 



G'p = [' dxK^ixh^olGciy + lnx)GF{y + ln(l - x)) - G^(y)], (6) 
Jo 



^In the forthcoming issues of the Physics-Uspekhi journal, it is intended to pubhsh the surveys on the effective action, 
on description of soft processes in quantum chromodynamics (physics of smah x and its relation to the Regge approach). 



where G'{y) = dG/dy, y = ln(p9/Qo) = ln(2Q/Qo) is the evolution parameter, defining the energy 
scale, p is the initial momentum, is the angle of the divergence of the jet (jet opening angle), assumed 
here to be fixed, Q is the jet virtuality, Qq = const , is the number of active fiavours. 



2 2iV,a5 



= (7) 

77 



the running coupling constant in the two-loop approximation is 



c^s{y) = i^{l-^-^-^]+0{y-'), (8) 




where 



the labels G and F correspond to gluons and quarks, and the kernels of the equations are 

K^{x) = ^ - {1 - x)[2 - x{l - x)], (10) 

KEix) = ^[x' + {l-xn (11) 



1 + - 
2 



(12) 



Nc=3 is the number of colours, and Cp = {N^ — l)/2Nc = 4/3 in QCD. The asymmetric form (|TOD 
of the three-gluon vertex can be used due to the symmetry properties of the whole expression. If one 
puts = in equation (||) and omits equation (|^), then one gets the equation of gluodynamics briefiy 
described above. The auxiliary variable u has been omitted in the generating functions for simplicity. 

The typical feature of any field theory with a dimensionless coupling constant (quantum chromo- 
dynamics, in particular) is the presence of the singular terms at x — in the kernels ([lOD , ([T^) of the 
equations. They imply the uneven sharing of energy between newly created jets and play an impor- 
tant role in the jet evolution giving rise to its more intensive development compared with the equal 
proportions (nonsingular) case. 

Let us note that these equations can be transformed into the linear equations for the moments of 
multiplicity distributions (see subsections 3.1 and 3.3). For the running coupling strength (H), they 
have been solved only in the perturbative theory approach. The systematical method, proposed in P0| , 
consists in using the Taylor series expansion of the expressions under the integral sign. It results in 
the modified perturbative expansion of the physically measurable quantities. Namely such solutions are 
usually considered when the analytical formulas are compared with experimental results. 

At the same time, these equations can be exactly solved ^ if the coupling strength is assumed 
fixed, independent of y, i.e. a5=const instead of (P). This becomes possible because of the scaling 
property leading to the power-like behaviour of the average multiplicity and all higher moments for the 
fixed coupling case. This behaviour differs from the dependences obtained in the case of the running 
coupling strength (see subsection 3.1). Therefore, the slow logarithmical decrease of this strength with 
transferred momenta increasing is crucial for the development of the parton cascades. 

Even though the system of equations (H), (P) is physically appeahng, it is not absolutely exact, i.e., 
it is not derived from the first principles of quantum chromodynamics. One immediately notices this 
since, for example, in these equations there is no four-gluon interaction term which is contained in the 
lagrangian of QCD. This interaction should, probably, correspond to the contribution to the integral 
term of the equations proportional to the product of three generating functions with a corresponding 
weight factor. Such a term does not contribute the singularity to the kernels and its omission is justified 



in the lowest order approximations. There were no attempts to take this term into account in the higher 
order perturbative expressions. 

Nevertheless, the modified series of the perturbation theory with three-parton vertices is well re- 
produced by such equations up to the terms including two-loop and three-loop corrections. As shown 
in Ref. 0, the neglected terms would contribute at the level of the product of, at least, five gener- 
ating functions. The physical interpretation of the corresponding Feynman graphs would lead to the 
treatment of the 'colour polarizability' of jets. 

These equations are justified up to some approximation of the modified perturbation theory (see 
below), because they only include those Feynman diagrams, where the gluons have strongly ordered 
transverse momenta, and are not close to the kinematical limit. Apart from this, they take into account 
the non-perturbative effects (e.g., the properties of the QCD vacuum) only in a simplified manner by 
the direct cut-off of the cascade evolution at some virtuality Qq. In principle, the effective infrared-safe 
coupling constant (without the Landau pole) may be used as the substitute for the phenomenological 
parameter Qq. It must be universal for different processes and tend to a constant limit at low virtualities. 



The constant average value of the coupling strength in this region has been used in papers |^ 
dealing with the non-perturbative corrections. However, the behaviour of the coupling strength is not 
the only non-perturbative effect. Therefore, we will use the more traditional perturbative approach. 
The non-perturbative effects in the three-jet events have been studied in more detail in the papers |^ . 



There are some problems also with the definition of the evolution parameter, with preasymptotic 
corrections etc. For example, the lower and upper limits of integration over x in Eqns (^), (|^) are 
constant and correspond to their asymptotic values. In reality, they vary in the preasymptotic region. 
Their form is determined by the restriction imposed on the transverse momentum which is given by the 
inequality 

kt = x(l - x)pe' > go/2. (13) 

This condition originates from the requirement that the formation time of a gluon {tfor-m ~ k/kf) must 
be less than its hadronisation time {thad ~ kR^ ~ k/Q^). It should be imposed for the perturbative 
QCD to be applicable. This leads to the requirement that the arguments of the generating functions in 
Eqns (I), (^ should be positive. Therefore, we must integrate in Eqns (|]), (^ over x from exp(— ?/) to 
1 — exp(— y). However these limits tend to and 1 at high energies [y — > oo). The omitted contributions 
decrease in the power-like manner with energy increasing. That is why it seems reasonable to learn 
more about the solutions of equations (|^), (^) near the asymptotic region taking first into account the 
perturbative (logarithmically decreasing) corrections, and only then take the neglected power-like terms 
into account as further corrections to these solutions. 

Moreover, this cut-off of the limits of integration is of physical importance. With the limits equal to 
exp(— y) and 1 — exp(— ?/), the partonic cascade terminates at the perturbative level Qo/2 as is seen from 
the arguments of the generating functions in the integrals. With the limits equal to and 1, one extends 
the cascade into the non-perturbative region with low virtualities Qi ~ xpQ/2 and Q2 ~ (1 ~ x)pQ/2 
less than Qo/2. Namely this region contributes terms of the order of exp(— y), power-suppressed in 
energy. It is not clear whether the equations and LPHD hypothesis are valid down to some Qq only or 
the non-perturbative region can be included as well. 

Some approximations are used to solve these equations with the running coupling strength. The 
Taylor series expansion leads immediately, as will be shown below, to the perturbation theory 
series in the exponent of the physically measurable quantities. It implies the summation of some specific 
subseries of the purely perturbative expansion in the coupling strength for this particular characteristics. 
This justifies the term of the "modified perturbative expansion" attached to this procedure. Moreover, 
the expansion parameter is not as itself but its square root 70. The asymptotic results (for extremely 
high energies) are obtained in the so-called double-logarithmic (DLA or DLLA if it is called as double- 
leading- logarithmic) or leading order (LO) approximation where the terms {as In^ s)" are summed. Here 
s is the cms energy squared. The emitted gluons are assumed to be so soft that the energy-momentum 



The corrections accounting for conservation laws in the G^G term, i.e., the shift in their arguments 
in equations (^, (|^), as well as the higher order terms in the weight asK (in particular, the non-singular 
terms of the kernels K and the dependence of as on the transverse momentum kt (|l5)) appear in the 
next-to-leading (NLO or MLLA - modified leading logarithmic, NLLA - next-to-leading logarithmic 
approximation) and higher (2NL0, 3NL0, ...) orders. Formally, these equations have been proven 
only for the next-to- leading (NLO) order of the modified perturbative QCD. However, one can try to 
consider them as the kinetic equations and solve them in higher orders with a hope to get from the 
obtained solutions and their comparison with experiment any indications on the role of the adopted 
assumptions. 

To have a guide to the uncertainties at higher orders of perturbation theory, one can try to generalize 
these equations by including the abovementioned effects in a more rigorous way than it is usually implied. 
However, there is no unique way of doing it. Also, it is not clear how one may modify them to include the 
non-perturbative effects, the colour coherence, the non-linearity at high densities etc, even though the 
very preliminary phenomenological versions have been proposed 48, 32, 49, 50, 51, 52, 33, 53|. The 



simplest proposal would be to compare two alternative evolution equations which use somewhat different 
assumptions leading to different higher order contributions as has been demonstrated in |^, |53 |. 

At the end of the section, let us mention that no direct solution of equations (|^), (|^) for the generating 
functions depending on both variables u and y has been obtained. Such solutions would be of interest 
for finding the location and the strength of the singularities of G in the complex plane u. This is 
important in connection with the behaviour of the moments of the distribution (see formulas ([26|), (pT]) 
below) and with some analogies from the statistical physics, where these singularities would indicate 
the phase transition point The corresponding singularities have been found only in the lowest 

order of the modified perturbation theory. It occured that they are located at the point u = 1 + Zq, 
where zq = C/{n){C ~ 2.552) tends to at high energies. The singular part of G at this location looks 
like 

G{u, y) = . f° + \ log ^^±1^ + 0(1). (14) 

[u — \ — Zo) u — 1 — Zq 3 Zq 

The structure of singularities is rather complicated. They are close to the point u = 1 where all 
moments are calculated. Moreover, with the energy increase they move closer to this point and coincide 
with it in asymptotics. There is, however, no special reason to worry about it since the generating 
function should be equal to 1 at this point according to its definition, and, therefore, all singularities 
must cancel somehow there. This is seen from the above expression as well. Nonetheless, this indicates 
the completely different structure of the singularities which one obtains in the lowest order perturbative 
expressions compared with the final result. 

In experiment, one always has to cut off the sum over the multiplicity in the definition of the 
generating function (|^) at some maximum multiplicity n^ax-, defined either by the energy- momentum 
conservation laws or by the definite conditions of a particular experiment. Then the generating function 
becomes a polynomial of the order of rimax with positive coefficients and, therefore, possesses just this 
number of the complex conjugated roots. Thus, one can find out the location of the singularity from 
experimental data only by increasing energy (and, consequently, the value of Umax) and following the 
evolution of their locations. According to the papers the complex-conjugated roots of the 

partition function of the grand canonical ensemble pinch the real axis of u at the limiting transition 
to the infinite volume just at the singularity location. The behaviour of the same type has been 
demonstrated in multiparticle production experiments for increasing rimax (in niore detail, see Ref. M). 



3 Comparison with experiment 

Let us turn now directly to the comparison of the theoretical results obtained with available experimental 
data. The main bulk of the data is provided by e+e^-processes at the energy. Many results can be 



found in the compilation [p^ ]. 

In experiment, one registers the particles created at the final stage, mostly, hadrons. Then the crude 
data are corrected for the effectiveness of the detectors, possible radiation of photons before the collision 
etc. At the same time, as has been mentioned before, the theory provides us with knowledge of the 
properties of quark and gluon jets and of the distributions of partons created during their evolution. 
Therefore, the problem of the correspondence of the theoretical results to the experimental data arises. 
One must separate the jets in an adequate manner and compile the "dictionary" for the translation 
from the parton to hadron " language" . 

In e^e~-annihilation processes, the quark-antiquark pair is first created. Therefore, the measurement 
of the multiplicity of the quark jet is of no problem because it is just twice smaller than the total 
multiplicity and equals to the multiplicity in one of the hemispheres. In terms of the generating 
functions this can be expressed by the relation 

Ge+e- ~ Gl. (15) 

To get the analogous results about the gluon jets, one should, however, have the access to the pairs of 
the gluon jets created by a colour-neutral source. This is necessary for the complete correspondence to 
the theoretical definitions. Unfortunately, such a separation of these events is possible with satisfactorily 
high precision only either in the decay processes T 'ygg 7+hadrons [0 or in the rare collisions 



leading to the almost parallel heavy quark and antiquark in the same hemisphere with the gluon in 
the opposite hemisphere |5^, ^ The high experimental statistics at Z" energy allows to do this. 
At other energies, where the number of registered events is much smaller, the separation of gluon jets 
with the help of methods, which do not depend on the chosen algorithm (the so-called "unbiased jets"), 
became possible quite recently after application of the special analysis of the two- and three-jet events 



according to two variables - the transverse momentum and energy - proposed in |62[ and used 



in (see the subsection 3.6). The separation of the quark and gluon jets is often done with the 



help of some special algorithms^. Their use is sometimes not completely identical to the theoretical 
requirements (the so-called "biased jets"). Then one has to rely only on the comparison with the results 
of those Monte Carlo models where the applied algorithm has been taken into account. The choice of 
the algorithm is by itself determined by its physical reliability, by the convenience for the analytical 
estimates and by the role of the hadronization corrections in these estimates. 

What concerns the hadronization stage, i.e., the transition from partons to hadrons, the Monte Carlo 
models use various phenomenological approaches mentioned briefiy above. In the analytical calculations, 
the variation of the parameter Qq in the combination with the hypothesis of the local parton-hadron 
duality usually plays the role of the "dictionary" and leads to reasonable results. 

3.1 The energy dependence of mean multiplicity 

The equations for the average multiplicities in jets are obtained from the system of equations (|^), (^ 
by expanding the generating functions in the power series of m — 1 and keeping only the terms with q=0 
and 1 with account of the definition of the average multiplicity (n) as 



u=l 



J2nP^ = {n). (16) 



dG 
du 

Finally, one gets the linear integro-differential equations with the shifted arguments under the integral 
sign for the average multiplicities. They read 



{nciy))' = I dx-f^[KS{x){{nG{y + Inx)) + {ndy + ln(l - x)) - {ndy))) 

+nfK^{x){{nF{y + Inx)) + (n^(y + ln(l - x)) - (ricd/)))], (17) 



'The most popular one is the method of the jet separation according to their relative transverse momentum called kt- 



dx^lKp{x){{nG{y + In a;)) + {npiy + ln(l - x)) - {np{y))'). 



(1^ 



From here, by solving these equations, one can learn about the energy evolution of the ratio of 
average multiplicities between gluon and quark jets r and of the QCD anomalous dimension 7 (the 
slope of the logarithm of average multiplicity in the gluon jet) defined as 



1 



{no)' 



They have been represented by the perturbative expansion at large energy (or large y) as 

7 = 7o(l - ai7o - ^270 - «37o) + ^(70), 

r = ro(l - ri7o - r27o " ^37o) + 0(70). 
Using the Taylor series expansion |^ of (n) at large y in Eqns (^), (|n 



(19) 

(20) 

(21) 



with (pOD, (pl|) and equating 
the terms of the same order in 70 in both sides, one gets the coefficients a,, Vi shown in the Table 1 
(ro = 9/4). 
Table 1 
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0.185 


0.426 
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0.280 


- 0.379 


0.209 
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0.191 


0.468 


0.080 


0.297 


- 0.339 


0.162 
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0.198 


0.510 


-0.041 


0.314 


- 0.301 


0.112 



The parameter 7 determines the exponent o 



the mean multiplicity of the gluon jet 



exp( 



i{y')dy'), 



(22) 



and therefore its perturbative expansion in terms of 70 corresponds in each particular order in 70 to 
the summation of some subseries in the purely perturbative expression for the average multiplicity, i.e., 
to the modified perturbation theory for this particular physical quantity. The choice of the gluon jet 
is related, first of all, with the tradition to study in the beginning the limiting case of gluodynamics 
where Uf = 0, i.e., quarks and, correspondingly, equation ([TB|) are not cosidered. This is justified also 
because in the lowest order approximations (the NLO corrections including) the energy dependence of 
average multiplicities in quark and gluon jets do not differ (for more detail, see next subsection). 

One of the most spectacular predictions of QCD states that in the leading order approximation 
(i.e., asymptotically), where 7 = 70, average multiplicities should increase with energy |^5|, |6^ 
like exp[2cVlog s] with the theoretically calculable value of c. This behaviour is just in between the 
power-like and logarithmical dependences correspondingly predicted by the hydrodynamical and multi- 
peripheral models^. Next-to-leading order results account for the term with ai in Eqn. (pO|) ||68| , |69| , [70 



and contribute the logarithmically decreasing factor to this behaviour. Namely these two terms deter- 
mine the main energy dependence of average multiplicities, and they are the same for quark and gluon 



jets. The higher order terms do not practically change this dependence |7T 

/3: 



72 



{no) = Ky-^^'^' exp{2c^ + —[2a2c' + ^{\n2y + 2)] 

2 '^1/^1 



r(ln% + i)])' 



(23) 



where c = (4A'"c//5o)^^^- The gluodynamics expressions can be obtained from this formula for = 



taking into account it also in the analytical expressions for given in 



^ It is not excluded, however, that the non-hnear corrections which appear, say, due to the high gluonic density 
can unitarize this approximation and give rise to the logarithmic increase of the mean multiplicity in the asymptopia 

rh 1I lorH lorJl 



The fitted parameters in the final expression are the overall constant normalization factor K which 
is defined by the confinement^ and the scale parameter Qq = 2A. The e"'"e~-data in the energy interval 
from the T-resonance to LEP-2 (i.e., approximately from 10 to 200 GeV) are well fitted by such an 
expression as seen in Fig. 1. The dotted line corresponds to the fit with two adjustable parameters. The 
dashed line shows the fit by the Monte Carlo model HERWIG. Other lines determined by the difference 
between gluon and quark jets will be explained in next subsection. 

Let us note here that the expansion parameter 7 is rather large at present energies ranging from 0.4 
to 0.5. The obtained expressions are valid at high energies and do not pretend to describe the energy 
behaviour of mean multiplicities near threshold. The computer solutions |j73|, |7^ of equations (|^), (H) 



give rise to the satisfactorily precise agreement with experimental data even at rather low energies. 

Equations (p!7D, (|T8|) can be solved exactly in the case of the fixed coupling constant because 
of the scaling property according to which the relation 

{niy + \nx)) ^ ^^4) 



Hy)) 

is valid for 7o=const and integro-differential equations (0) are reduced to the system of easily 

solvable linear algebraic equations. The energy increase of mean multiplicity becomes power-like because 
the exponent acquires the general structure of the type of 70 Ins. 

3.2 Difference between quark and gluon jets 

The system of two equations for quark and gluon jets predicts that asymptotically the energy dependence 
of mean multiplicities for them should be identical. Moreover, this coincidence is exact in the next-to- 
leading approximation of the modified perturbation theory as well. Higher order correction, even though 
violating this beautiful feature, are comparatively weak in the functional dependence on energy. Namely 
this explains the initial success in the description of the energy dependence of average multiplicities in 
e"'"e~-annihilation in the framework of gluodynamics. 

The absolute normalization is not fixed. However, the relative normalization, as given by r is 
calculable. The gluon jets are more "active" than the quark jets so that the ratio r = (jig) / {np) of 



average multiplicities between gluon and quark jets should tend at high energies |]75| to the ratio of the 
Casimir operators Ca/Cf = 9/4. 

Once again, by comparison with experiment this prediction shows how far are we now from the true 
asymptotics. Even though this prediction obtained in the leading order (LO) approximation is valid 
qualitatively, its quantitative value is still rather far from experimental ones where this ratio is about 
1.5 at energy and even smaller at lower energies. The higher ordeiQ terms |76|, |78|, |6^, |72|, |79 



are rather important just for this parameter. They have been calculated now up to 3NL0 (or 4NL0; 
see the footnote) terms (see Table 1) and improve the agreement approaching the experimental value 
with an accuracy about 15% (see Fig. 2 where the old notation for the approximation order has been 
used). Let us stress again that just the ratio of mean multiplicities r (but not the energy behaviour of 
any of them, {no) or (n^), separately) is most sensitive to these corrections. This is because the main 
LO and NLO terms cancel there and decline from the constant value tq is governed by the higher order 
terms. This decline can provide a guide to further generalizations of the equations and to the proper 
account of the nonperturbative contributions. 

Surely, the higher order terms change slightly also the energy behaviour of multiplicities for quark 



^That is why the lower hmit of the integration over y' in formula ( p2[) is not fixed. 

^The difference in the definition of the "order" for the anomalous dimension 7 and for the ratio of mean multiplicities 
r is described in Ref. f!% . In particular, the term in r should be actually considered as 4NL0 correction since it is 
added to the still uncalculated term 04 in the higher order expression for the quark jet anomalous dimension and, 
consequently, influences its multiplicity. 



jets compared to gluon jets[| as observed in experiment. However, the simultaneous fit of quark and 
gluon jets witli tlie same set of fitted parameters even in tfie framework of 3NL0 approximation is still 
not very accurate as is seen from the shaded area in Fig. 1. This area demonstrates the attempt of such 
a description with the data on gluon jets divided by the theoretically (analytically) calculated value of 
the ratio r. Its failure shown by the shift of the shaded area compared to the experimental points and 
by its large width is again due to the insufficiently precise description of the ratio r. The agreement 
with experimental multiplicity in e+e~ is restored either in the case of the normalization according to 
gluon jets with the only adjusted parameter Qq = 2A (the solid line in Fig. 1) or in the case with both 
parameters fitted (the dotted line). 

The computer solution of the equations gives rise to very good agreement with experiment at 



the Z° resonance and leaves a rather small difference about 20% even at so low energies as the mass 
of the T. In the analytical approach, the quite good agreement on the ratio r at the energy of the 
has been achieved when the equations are modified to account for the phase space limitations imposed 



by energy-momentum conservation in the dipole cascade picture of the string approach |52|. However, 
some problems arise for higher moments of the multiplicity distribution in such an approach ||53[] . 

Equation (|2T| ) states that the ratio r increases with energy increasing due to the decrease of 70 (the 
asymptotic freedom) but its slope becomes smaller. Experimental values also increase with energy as 
seen from Fig. 2. There are common tendencies and good qualitative agreement between the results of 
analytical calculations and experiment on the value of the ratio r. However, the quantitative statements 
are valid only within the 15 - 25% accuracy. 

More detailed characteristics of jets behaviour can be obtained from the slope and the curvature 



of energy dependences of multiplicities [^, |72[. Some experimental data about it have been recently 
obtained ||63l . 

The slope of the energy dependence of the ratio of average multiplicities between gluon and quark 
jets r, i.e., the derivative r' ("the slope of the ratio"), decreases with the energy increase according to 
the expansion (^ as ln"^/2 s. It is very sensitive to the higher order corrections |^D|, |7^, first of all, 
due to the large value of the ratio r2/ri: 

r' = Br.ralll + ^ + + ^iH + 0(7o')], (25) 

where B = Pq/SN^, Bi = Pi/ANcPo- Therefore its use for comparison with experiment does not seem 
to be informative enough. 

At the same time, such corrections are partly cancelled in "the ratio of the slopes" r*-^-* = (^g)V (^f)' 



for gluon and quark jets |80, Thus, the ratio of the slopes must depend on energy much weaker. 



This has been observed in experiment. The first comparison P] of the theoretical estimates of this 



ratio I^T], ^] with experiment has shown quite good agreement, however, still within rather large 
undefiniteness in experimental results. Recent more accurate separation of unbiased gluon jets 
made it possible to get the values not only the ratios r and r^^-* at energies from T to Z° (including the 
intermediate ones) but also the values of the ratio of the second derivatives r^^^ = (nc)" / {np)" . The 
asymptotic values for all of them are predicted equal 2.25. The corresponding values at the energy of the 
hadronic component in ^^-system 30 GeV are r = 1.422 ± 0.051, r^^) = 1.761 ± 0.071, r^^) = 1.92 ± 0.13 
and at 80 GeV r = 1.548 ± 0.042, r(^) = 1.834 ± 0.09, r^^) = 2.04 ± 0.14. The error bars are mainly 
determined by the systematic errors. These results are in agreement with the analytical QCD predictions 
5D|, |72[ that at present energies r^^^^ should be closer to its asymptotic value of 2.25 than r^^\ and r^^^ 



closer to 2.25 than r, i.e. the inequalities r < r^^^ < r^^^ are valid. The experimentally obtained values 
of the ratios of slopes and curvatures of energy dependence of gluon and quark jets and their comparison 
with some theoretical calculations are demonstrated in Figs. 3 and 4. 



^The value 7 in (22) is replaced there by 7^^ — 1^ f' /fi which differs from 7 only in higher orders as is easily estimated 
if one recalls that r' ^ "f^ ^ 7q (see below formula (|2F 



One can hope that the asymptotic values of the multiphcities ratio r can be approached closer if 
the soft particles are only considered with energies much lower than the total energy of the jet |84|. 
This effect has been really noticed but it is not very strong. The ratio r for particles with low momenta 



is about 1.8 85 



It is interesting that in the Monte Carlo model HERWIG accurately reproducing this ratio for soft 
hadrons the corresponding ratio for soft partons occurs exactly equal to 2.25. Thus the role of the 
hadronization stage is very important for the ratio r in this particular model. Let us note that this 
situation differs from the case where all momenta are averaged as demonstrated in Fig. 2 with the lines 
of the HERWIG model practically indistinguishable for partons and particles. This is related to the 
strong difference of the spectra which are rather soft for hadrons and much harder for partons in this 
model (for more detail, see Ref. |^). The direct use of the local parton-hadron duality is, therefore, 
impossible here because of the large value of the adopted cut-off parameter Qq indicating the early 
finish of the parton stag^ 

The parton stage in the evolution of the jet of final hadrons can be separated in a following way. 
In the HERWIG model it is assumed to finish at Qo=0.8 GeV. Thus it corresponds to the rather large 
transverse momenta. Therefore one should also choose the hadrons with large transverse momenta. 
Analysis of the data about the comparatively soft hadrons (p < 4 /c) with rather large transverse 
momenta pt >0.8 /c have led to the result |^5[ for the ratio r = 2.29 ± 0.017 in this region which 
agrees with asymptotic predictions. In its turn, this result could agree with the computer solutions 
of equations (||), (H) and with the local parton-hadron duality hypothesis only if one adopts that the 
evolution of the "partonic" cascade in the computer calculations to the region Qq < 0.8 GeV plays the 
same role as the "hadronization" stage in the HERWIG model, where such an evolution is discarded. 

The widths of the multiplicity distributions differ for quark and gluon jets, the former being some- 
what wider. Qualitatively, QCD describes this tendency but quantitative estimates are rather uncertain 
yet as is discussed in more detail in subsection 3.8. 



3.3 Oscillations of cumulant moments 



The shape of the multiplicity distribution can be described by its higher moments related to the width, 
the skewness, the kurtosis etc. Equations for the moments of the distributions of the g-th rank are 
obtained from the system of equations (||), (H) by comparing the terms of the type {u — 1)'^ in both sides 
of equations. The g-th derivative of the generating function corresponds to the factorial moment Fg, 
and the derivative of its logarithm defines the so-called cumulant moment Kg. The latter ones describe 
the genuine correlations in the system irreducible to the products of lower order correlations (this recalls 
the connected Feynman graphs). 



P„n(n — l)...(n — q + 1) 
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These moments are not independent. They are connected by definite relations which can easily be 
derived from their definitions in terms of the generating function: 



q-l 



Fq ^ ^ Cg_iKq—mFiji. 
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m=0 



They are nothing other than the relations between the derivatives of a function and of its logarithm at 
the point where the function itself equals 1. Here 

(g-1)! r(g) 1 
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q-l 



m!(g — m — 1)! r(m + l)r(g 
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mB{q, m) 



(29) 



are the binomial coefficients, and F and B denote the gamma- and beta-functions, correspondingly. 
Thus there are only numerical coefficients in recurrence relations (^) and the iterative solution (well- 
suited for computer calculations) reproduces all cumulants if the factorial moments are given, and vice 
versa. In that sense, cumulants and factorial moments are equally suitable for analysis of multiplicity 
distributions. The relations for the low ranks are 

Fi = Ki = l, 
F2 = K2 + I, 

F3 = KS + 3K2 + I. (30) 

The Taylor series expansion of equations (^, (^ gives rise to the term G{y) in the right-hand 
side. Dividing by it the left-hand side G'{y) one gets the derivative over y of the logarithm of the 
generating function. The gluodynamics equation in the leading order approximation can be reduced to 
the differential equation of the second order: 

(lnG(y))" = 7o'(G'(i/)-l); G{0) = u, G'(0) = 0. (31) 

After that, the derivatives over u of both sides considered at u = 1 lead in a natural way to the prediction 



40| , PI , |79| concerning the behaviour of the ratio Hq = Kq/Fq. At asymptotically high energies, this 
ratio as a function of the rank q is predicted to behave^ as q~^. Thus the role of genuine correlations 
is strongly damped in systems with the large number of particles. 

However, the asymptotics is very far from our realm. At present energies, according to QCD, this 
ratio should reveal the minimum at g ~ 5 and subsequent oscillations. In gluodynamics, the minimum 
is located at 

24 

to ~ T— + 0.5 + 0(70). (32) 
II70 



This astonishing qualitative prediction |^ of the new type of the moments behaviour has been 
conffimed in experiment (for the very ffist time in Ref. [Q) as in e~^e~ (see Fig. 5 [^) as in 



hadronic processes PSI . The predicted negative minimum of Hq is clearly observed. These oscillations 
can correspond to the replacement of attractive forces (clustering) by repulsion (between clusters) in 
systems with different number of particles. Let us note that the colour coherence in QCD discussed 
below also predicts both the attractive and repulsive forces in the systems of the coloured partons. 

With energy increasing, one should observe ffist the disappearence of the oscillations at high ranks 
and then the smaller dip at the first minimum and its slow shift to higher values of q from its initial 
position at g ~ 5 at Z° energy due to the decrease of 70 (see Eqn ([3^)). Finally, the ratio Hq will slowly 
tend to its asymptotic dependence 

The minimum position slowly changes with energy and with the size of the phase space window 



because it is inverse proportional (see [^, |SD|) to the square root of the running coupling strength, i.e., 
to 7o. For some specific processes this shift can be very strong (e.g., it has been found for instanton 
induced processes that the minimum moves to g ~ 2 because the multiplicity distribution in these 



processes is very narrow). At the same time, other experimental data conffim that the oscillations 



reveal themselves on the qualitative level in various processes, i.e., have some universal origin |^8|. The 
graphical presentation of the ratio Hq in place of the moments is suitable, because the moments increase 
fast with the rank q increasing but this growth cancels in their ratio. 

The quantitative analytical estimates are not enough accurate because, ffist of all, the expansion 
parameter becomes equal to the product 57 which is close to one or even exceeds it for all g > 1. 
Therefore the perturbative approach is, strictly speaking, inapplicable to this problem. However, some 
tricks like Fade-approximation can be used to improve it. At the same time, the numerical computer 
solution |73, WM reproduces oscillations quite well. 



This is clearly seen from equation ( 31 ) because each differentiation over y in the left-hand side provides the factor 57 



These new dependences differ from all those typical for the previously attempted distributions of the 
probability theory, in particular, such previously popular in particle physics distributions as the Poisson 
and negative binomial ones. The characteristics which are special for the newly found distribution at 
high energies have been discussed in Ref. 



31 



3.4 The hump-backed plateau 



Dealing with inclusive distributions of partons within the jet over the rapidity or over the share of the 
jet energy (x) acquired by the parton, one should solve equations of the type for the generating 
functional. The variation over the auxiliary function u{k) gives rise to integro-differential equations for 
the one-particle distributions and the correlation functions. In the leading order in gluodynamics the 
equation for the distribution function D{x,y) is written as 



d{xD{x,y)) 
dy 
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(33) 



with the boundary condition xD{x,0) = S{1 — x). Here, x = k/Ej, k is the particle momentum, Ej 
is the jet energy, K{z) = Kq{z). In QCD there are two equations for gluon and quark jets. The 
subsequent use of the Mellin transformation allows one to reduce these equations in the low orders to 
the solvable differential equations. This has been done up to the NLO approximation (see 0, ||). 

As predicted by QCD, the momentum (or, more accurately, rapidity yJQ) spectra of the particles 
inside jets in e+e~-annihilation processes should have the shape of the hump-backed plateau 
|67| , pO[| . This striking prediction of the perturbative QCD differs from the previously popular flat 
plateau advocated by Feynman. This has been found in experiment (Fig. 6). The depletion between 
the two humps is due to the angular ordering and colour coherence in QCD. The humps are of the 
approximately Gaussian shape (the distorted Gaussian) near their maxima if the variable ^ = In - is 
used. This prediction was first obtained in the LO QCD, and more accurate expressions were derived 
in NLO 0. 

The maximum position in the ^-distribution must shift with energy increase according to the 
almost logarithmical law with some corrections of the type 0(ln~^''^ s): 



^o = y 
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Moments of the distributions up to the fourth rank have been calculated. The spectrum drops quite 
fast towards smaller momenta. This is a consequence of the colour coherence for the ever softer gluons. 
It becomes especially noticeable if the variable ^ is used. In the absence of the colour coherence the 
shape of the distribution near the maxima losts its Gaussian form, the maxima are shifted to the region 
of smaller values of x (larger ^) and, therefore, the spectrum becomes much softer. 

The comparison of these QCD predictions with experimental data at different energies has revealed 
good agreement both on the shape of the spectrum (see Fig. 7 for e^e~ from IQ) and on the energy 
dependences of its peak position (see Fig. 8 for e~^e~, ep,pp from |^) and of its width. 

For soft particles, the spectrum does not depend on the energy and has the shape p3|, p^: 

, 1? , « cxsih). (35) 
dyed m kt 

This explains the discussed above possibility to observe the approach to the asymptotic value of the 
ratio r even at finite energies if the soft particles with pt —>■ Qo are chosen. The physical origin of this 
effect lies in the ability of soft partons to resolve only the total charge of the jet (i.e., of the parton 
initiating it) but not its internal structure. 



3.5 Difference between heavy- and light-quark jets 



Another spectacular prediction of QCD is the difference between the spectra and multiphcities in jets ini- 
tiated by heavy and hght quarks. Quahtatively, this corresponds to the difference in the bremsstrahlung 
by muons and electrons where the photon emission at small angles is strongly suppressed for muons 
because of the large mass in the muon propagator. Therefore, the intensity of the radiation is lower in 
the proportion to the ratio of the masses squared. The suppression of the photon radiation at small 
angles was proposed to be exploited also for the top quark where it is especially strong due to its 
large mass. However, the intensity of the radiation would be too low due to the same reason. 

The coherence of soft gluons also plays an important role in QCD. For heavy quarks the accompa- 
nying radiation of gluons should be stronger depleted in the forward direction (the so-called dead-cone 
or ring-like emission). At large angles there should be no visible difference between the emission of 
gluons by heavy or light quarks. In spite of the close analogy, there is a strong difference between 
quantum chromodynamics and quantum electrodynamics connected to the fact that the gluons possess 
the self-interaction while the photons have no point-like interaction of one with another. 

It was predicted |100| , |101|] that the coherence of the soft gluons should result in the energy- 



independent difference of companion mean multiplicities for heavy- and light-quark jets of equal energy. 
The companion hadron multiplicity for the heavy quark is defined as the difference between the total 
multiplicity of the jet initiated by this quark and the constant decay multiplicity n'^'^ in its rest system 
(which is equal to 11.0±0.2 for the 6-quark and 5.2±0.3 for the c-quark). Correspondingly, the difference 
of the total average multiplicities n is written as 

nQQ{E) - n,,{E) = n^Q{M) - n,,{^eM)[l + 0(«s(M))]. (36) 

Here, g, Q denote the light and heavy quarks, E are their energies, M is the heavy quark mass, \/e ~ 
1.65. The naive model of energy rescaling ||102| , |103|| predicts the decreasing difference at higher energies. 
The experimental data (see Fig. 9 from ||104]| ) comparing the multiplicities in jets initiated by light 
(u, d, s) and heavy (6, c) quarks support this QCD conclusion. 

Another consequence of such a distinction is the effect of the leading heavy quarks, i.e., the con- 
sumption by them of the rather large share of the primary momentum in the process of the emission 
of gluons. This is related to the suppression of the emission of hard gluons with low transverse mo- 
menta which results in the low energy loss by heavy quarks. This effect has also been observed ||105|| in 
e"'"e~-annihilation processes. The inclusive spectra of the 6-quark jet at the energy of the have the 
maxima at xq ~ 0.8 — 0.9. In the hadronic processes, this effect was earlier used ||106|| to explain the 
so-called "long-flying cascades" in cosmic rays at the energy about 100 TeV. In result, the conclusion 
about the strong increase of the total cross section for the production of heavy quarks in the energy 
interval up to 100 TeV was obtained. It finds now the direct support from accelerator data. 

Also, it has been noticed ||107|| that the momentum spectra of particles in the 6-jets are much softer 
than those in the jets initiated by light quarks as has been predicted in QCD ||9^ with account of the 
coherence. 

Concerning the difference of the angular distributions for gluon emission by heavy and light quarks 
initiating the cascade, which lies at the background of all these effects, there exist just the preliminary 



results ||108| , |104|| , which favour the QCD predictions. 

Let us note that for hadronic processes the ring-like structure of the polar angle distribution can be 
observed also for the jets themselves, i.e., for the partons initiating these jets. It can result due to the 



limitations imposed on the radiation length or due to the so-called Cherenkov gluon effect ||109| , |110 
what corresponds in the "operationalists" language to the emergence of some effective mass. 



3.6 Colour coherence in 3-jet events 



When three or more partons are involved in a hard interaction, one should take into account colour- 



angles between the two jets |Q. For example, for two jets moving quite close to each other, one should 
consider their mutual screening. It implies that for the "resolution" of such a pair one should use the 
rather "hard probe" because the soft ones will feel only the total colour charge of the pair as a whole. 
In particular, it results in the "colour transparency" of such a pair during its motion inside a hadronic 
medium. In electrodynamics, the similar effect of the mutual screening of the electric charges of the 
closely moving electron and positron until they separate at a large distance is known as Chudakov effect 
ilTTI. 



In equation for the generating functionals one should take into account those correlations between 
the jets which arise if the Feynman diagrams are considered. However, the probabilistic interpretation 
in the framework of equations (|^), (|^) in some approximation may fail for these effects. Therefore, for 
studies of the colour coherence effects the Feynman diagram technique is preferred. 

Several such effects have been already observed. In particular, the total multiplicity can not be 
represented simply as a sum of flows from independent partons. In terms of the observed multiplicities 
rie+e- in annihilation processes, one can write down the multiplicity in the three-jet events nqqg as 

^^999 = K+e-i'^Eq) + 0.5r(pt)ne+e-(pt)](l + 0{as)), (37) 

where Eq is the quark energy and pt is the transverse momentum of the gluon in the center of mass 
system for the pair qq. The ratio of the multiplicities for gluon and quark jets is considered at the energy 
equal to the transverse momentum of the gluon which, in its turn, is related to the virtuality of the 
quark emitting this gluon. The existence of two scales in the three-jet events is clearly demonstrated in 
this way. The comparison of formula (0) with experimental data |p.l2| , |113|| has shown that it is valid 
for large angles between the jets. 

QCD predicts that the particle flows should be enlarged in the directions of emission of partons and 
suppressed in between them. Especially interesting is the prediction that due to the negative interference 
this suppression is stronger between the gg-pair than between gq and gq in the hard e"'"e~ — > qqg event 



if all angles between partons are large. This phenomenon is known as the "string" [[L14|| or "drag" 
115 | effect. All these predictions have been confirmed by experiment (see Fig. 10 from |p.l(j|| ). In qqg 
events the particle population values in the qg valleys are found larger than in the qq valley by a factor 
2.23±0.37 compared to the theoretical prediction of 2.4. Moreover, QCD predicts that this shape is 
energy-independent up to an overall normalization factor. 

Let us note that for the process e^e~ qq'j the emission of additional photons would be suppressed 
both in the direction of a primary photon and in the opposite one. In contrast, in the case of the emitted 
gluon one observes the string (drag) effect of enlarged multiplicity in its direction and stronger suppres- 
sion in the opposite one. This suppression is described by the ratio of the corresponding multiplicities 
in the qq region 

= (38) 

which is found to be equal 0.58±0.06 in experiment whereas the theoretical prediction is 0.61. 

The colour coherence reveals itself as inside jets as in the inter-jet regions. It should suppress 
both the total multiplicity of qqg events and the particle yield in the transverse to the qqg plane for 
decreasing opening angle between the low-energy jets. When the hard gluon becomes softer, colour 
coherence determines, e.g., the azimuthal correlations of two gluons in qqgg system. In particular, 
back-to-back configuration {ip ~ 180°) is suppressed by a factor ~ 0.785 in experiment, 0.8 in HERWIG 
Monte Carlo and 0.93 in analytical pQCD. 

Let us stress once again that the colour coherence determines the topological dependence of the 
jet properties predefined theoretically in terms of the parton diagrams. The interference between the 
gg-pairs not connected by colour is suppressed by the factor This interference is not accounted by 

the models using the Monte Carlo method. Thus their success in the description of experimental data 
implies the smallness of such effects. Nevertheless, these small colour-suppressed effects disappearing 



at A'^c ^ oo can become really important for the distinction between the analytic diagram approach 



in QCD and the purely probabilistic Monte Carlo schemes ||117|| . The study of the individual events 
with very high multiplicity may be crucial for getting the decisive conclusions. The most important 
lesson derived from the correspondence between the theory and experiment in this case consists in the 
conclusion that the colour coherence leads to the effects observed at the hadron level and is not wiped 
out by the hadronization stage. The hadron distributions depend on the topology of the parton stage 
with hard colour objects that gives further support to the local parton-hadron duality hypothesis. 

Some proposals have been promoted for the modification of formula ( P7D with more correct account 
for the phase space and for special two-scale (the energy and the transverse momentum) analysis of 
3-jet events when the restriction on the transverse momentum of a gluon jet is imposed 
This corresponds to the simultaneous account for the energy and the virtuality of the initially produced 
quark and antiquark. The first approbation of this proposal ||118|| has shown the correctness of such a 



modification. This is related to the solution of the problem of the unique separation of the three-jet 
events according to the kinematics at the parton level and clear detection of the gluon jet in this system 
mentioned above. In general, jets are biased because event selection according to some algorithms 
introduces bias on multiplicities. Jet properties depend on two scales, the transverse momentum and 
available rapidity range. To compare with theory, one should get unbiased results. For gluon jets, they 
have been obtained in using the following formula 



riggipt) = 2[nggg{s,pt) - ngg{s,pt)], (39) 
where the gluon transverse momentum is given by 

PtS = SqgSgg (40) 

and Sig is the squared cms energy of the system ig. This multiplicity depends on a single variable, 
whereas the terms in the right hand side depend also on the energies. This implies that there is no 
dependence on the algorithm of jets separation. It is this accounting of the dependence on the cutoff of 
the transverse momentum, below which the gluon jet is not resolved and the whole system is treated 
as a gg-pair, that finally leads to the independence of the multiplicity of the gluon jet on the adopted 
algorithm of jet separation. 

Let us note that the fast divergence of the modified perturbation series with increase of the term 



order ra, which is like n!, usually related to the notion of renormalons ||119|| , gives rise to large non- 
perturbative corrections. This divergence was noticed above as the appearance of 57 factor as the 
expansion parameter for high ranks of the moments of the distributions. 

The non-perturbative corrections to event shapes (especially, to three-jet events) should be men- 
tioned here. Just in this case these corrections are very strong, and there are numerous experimental data 
here. In distinction to the logarithmic dependences of the perturbative approach, the non-perturbative 
terms decrease in a power-like manner with increase of the transferred momentum. They can be written 



in a universal way [ 120 for various characteristics Vi as 

M(XQ(Ae//(p)/Q')^ (41) 

where the parameter Ae// does not depend on a particular characteristics of the event shape Vi, which 
can be the jet mass, thrust, jet widening etc. The numerical value of the parameter Cj is calculated in 
QCD [g. 

Even though such an approach is quite successful in more precise description of experimental data, 
there are many unsolved problems. In particular, it is not clear where one should cut off the pertur- 
bative series and add to it the non-perturbative terms, what is the relative role of corrections with 
increasing values of p at lower transferred momenta and where, in general, the non-perturbative 
terms (ETI) should saturate and flatten out. Moreover, it was supposed WT\\ that one should use the 



renormalization group improved perturbative expansion and it was shown that this improvement leads 
to very satisfactory description of experimental data. 

Results for jets in ep and pp processes also favour the theoretical expectations of the role of coherence 
effects for emission of soft gluons. Here, analysis is complicated by the internal structure of the colliding 
objects. Especially interesting are the data ||122| , |123|| on the topology of the pp-events either with the 
production of the IV-meson in combination with the quark or gluon jet or with the production of two 



hadronic jets |83, 124 1 



3.7 Intermittency and fractality 

The self-similar parton cascade leads to special multiparton correlations. Its structure with "jets inside 
jets inside jets..." caused by the angular ordering has provoked the analogy with turbulence and the 
ideas of intermittency [|125|| according to which the increase of fluctuations in ever smaller phase space 



volume studied (e.g., in smaller rapidity intervals 5y) must lead to the increase of factorial moments 
according to the following power-like behaviour: 

F, oc {6y)-'^'^''\ (42) 

In its turn, this self-similar structure should result in a definite geometric pattern. Namely, this leads 
to the fractal distribution of particles inside their available phase space |[126|| . The notion of fractality 



allows one to quantify the characteristics of the process expressing them in terms of the fractal (or 
multifractal, Renyi) dimensions. Especially important is the fact that the geometric dimensions are 
related to the physical characteristics, namely, to the intermittency exponent 0(g), which shows the 
slope of the increase of the factorial moments on the doubly logarithmic scale (see Eqn (0)). The 
slopes (the intermittency exponents) 0(g) for different ranks q are related to the Renyi dimensions Dg 
in a following way: 

0(g) = (g-l)(D-D,), (43) 

where D is the topological dimension of the analyzed phase space windows (for example, D=l if the 
dependence of the factorial moments on the length of the rapidity interval is studied). 

To calculate these characteristics, one uses the diagramatic approach in the same way as it was 
described above when discussing the colour coherence effects. This is necessary now because one has 
to deal with a small part of the total number of the created partons within the fixed small phase space 
volume. The fractal behaviour is usually defined by the dependence of the logarithm of the factorial 
moments in function of the logarithm of the size of the chosen phase space region. The (mono)fractal 
behaviour would display the linear dependence of logarithms of factorial moments on the logarithmic 
size of phase space windows. The moments are larger in smaller windows, i.e. the fluctuations increase 
in smaller bins in a selfsimilar power-like manner if the (mono)fractal distribution is studied (see the 
review paper [127]). 



In QCD, the power dependence of moments on the size of the phase space window appears for a 
fixed coupling regime In this case, the monofractal behaviour with a constant Renyi dimension is 
pronounced: 

D, = (44) 

and, correspondingly, the intermittency exponent 0(g) = const at any fixed rank g. 

The running coupling strength has a definite scale and, therefore, it leads to some decline from this 
simple self-similar (monofractal) behaviour. The running property of the coupling strength in QCD 
flattens off |p.28| , |129| , |13CI|| this dependence at smaller bins, i.e. the multifractal behaviour takes over 
there. 

Both the linear (in the doubly logarithmic scale) increase at comparatively large but decreasing 

Kins anH its flattpnincr for vprv small Kins havp> Kfpn nKsprvpH in p>Yr>primp>nt f spp> Pier 11 from flTSlH^ 



However, only qualitative agreement with analytical predictions can be claimed here. The higher 
order calculations are rather complicated and the results of LO with some NLO corrections are yet 
mostly available. In experiment, the different cuts have been used which hamper the direct comparison. 
However, those Monte Carlo models where these cuts can be done agree with experiment. The relative 
role of the partonic and hadronization stages of the cascade in this regime as well as the applicability 
of the local parton hadron duality hypothesis to the correlation characteristics of the process are still 
debatable. 

The transparent interpretation of the observed effect in terms of the fractal phase space volume has 



been proposed in the framework of the Lund dipole cascade model |132, 133 



3.8 The energy behaviour of higher moments of multipHcity distributions 

Differentiating both sides of equations (|), over u and using formulas (|26|), (|27|), one can get the 
equations for the moments of multiplicity distributions of any rank. Their solutions would describe 
the behaviour of the multiplicity fluctuations. They tell us that the fluctuations of the multiplicity 
of individual events must be larger for quark jets as compared to gluon jets. Therefore, the moments 
of their distributions in quark jets are larger than the corresponding moments for gluon jets. This 
tendency is clearly seen in experimental data ||116|| . The factorial moments increase both with their 



rank and with energy increasing. In the perturbative expansion one gets the formulas ||134|| , similar to 
formula (pT]) for r, and the energy increase of the moments is determined by the decrease of 70. The 
corresponding coefficients in front of the terms 79 are calculated in [[134|| . From the mentioned above 



behaviour of if^-moments one easily guesses that the same is true for the cumulant moments with some 
difference. 

The experimental results for the second rank factorial moments of 41.8 GeV gluon jets = 1.023 
and for 45.6 GeV uds quark jets, F2 = 1.082 are much smaller than the known long ago (see, e.g., [|l|) 
asymptotic predictions, viz. 1.33 and 1.75, respectively. The NLO terms improve the description of the 
data compared to the leading order results. If one accepts the effective value of averaged over all 
the energies of the partons during the jet evolution to be as ~ 0.2, one obtains [ p.34|| the NLO values 

^ 1.039 and ^ 1.068 at these energies which are quite close to the experimental results. In this 
sense the NLO prediction can be said to describe the widths of the gluon and quark jet multiplicity 
distributions at the energy to within 10% accuracy. 

Unfortunately, the 2NL0 and 3NL0 terms worsen the agreement with data compared to NLO (but 
not compared to LO) results ||134|| . The same is true for the higher order moments. It raises the general 



theoretical problem of the convergence of the perturbative expansion in view of the large expansion 
parameter qj mentioned above. The attempts to account for conservation laws more accurately by the 



modified evolution equations for high moments 53 1 have not led to the success yet. It is remarkable 



that the computer solution of the QCD equations provides a near-perfect description of the higher 
moments as well. This suggests that the failure of the analytical approach at higher orders is mainly a 
technical issue related to an inadequate treatment of soft gluons and of energy-momentum conservation. 



This conclusion is also supported by the success of the Monte Carlo model ARIADNE [ |135|| in describing 
these characteristics. 

The rather accurate experimental results about the multiplicity moments of the separated gluon jets 
136|| are available now up to the rank 5, and for the quark jets even up to the rank 17, because 



the accuracy of measuring the multiplicity distributions for the latter is much higher due to the larger 
statistics of the data. 

The higher moments of the multiplicity distributions are determined by the integrals of the corre- 
lation functions for partons inside the jet which depend on their angular distributions. The angular 
correlations have been studied as well ||137|| . The comparison with experimental data | |138| , p.39| | shows 
that the agreement becomes better with the energy increase even though it is still unsatisfactory at 
small angles. Some "infrared" stable characteristics like correlations of the energy and multiplicity flows 



have been considered (see, e.g., ||140| , |141|| ). They do not require for the determination of the jet axis, 
as in the case of the angular correlations, where the accuracy of this determination is crucial for the 
correlation studies at small angles. 

In the case of the fixed coupling constant, these equations are exactly solvable [H, |2|. Their 
solutions as functions of energy possess the scaling property. All moments behave in a power-like 
manner with energy increasing. 

Let us make two short technical comments at the end of this subsection. All analytical expressions 
should satisfy the requirements of the limiting case of the supersymmetric QCD (the mote detailed 
discussion of this problem and the corresponding formulas for the coefficients and rj from the Table 1 
see in Ref. |]T2|). The factorial moments are always positive according to their definition. These two 
requirements are the necessary but not sufficient conditions for the correctness of the calculations. 
However, the moments calculated in Ref. do not satisfy the second condition. At the same time. 



the coefficient r2 obtained in Ref. ||7^ satisfies the requirements of the supersymmetric QCD (as well 



as that from Ref. ||72[ which differs numerically) but the method of the renormalization group used in 
Ref. does not take into account correctly the energy-momentum conservation laws and gives rise 
to the numerically smaller value of this coefficient. 



3.9 Subjet multiplicities 

A single quark-antiquark pair is initially created in e"'"e~-annihilation. With very low angular resolution 
(large angle averaging) one observes two jets. A three-jet structure can be observed if a gluon with 
large transverse momentum is emitted by the quark or antiquark. However such a process is suppressed 
by an additional factor as, which is small for large transferred momenta. Its probability can be cal- 
culated perturbatively. At relatively low transferred momenta, the jet evolves to the angular ordered 
subjets ("jets inside jets inside jets..."). Different algorithms have been proposed to resolve subjets. 
By increasing the resolution, more and more subjets are observed. For very high resolution, the final 
hadrons are resolved. The resolution criteria are chosen to provide the infrared safe results (see, e.g., 
the papers H |4|, |T4|, pij |l45l |146|). 

In particular, one can predict the asymptotic ratio of subjet multiplicities in 3- and 2-jet events if 
one neglects the soft gluon coherence: 



SI 
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(45) 



Actually, the coherence reduces this value to be below 1.5 in experiment for all acceptable resolution 
parameters. The theoretical predictions agree quantitatively with experimental findings | |147| , |148| ] 
at the lowest resolution, where this ratio is equal to 3/2, and only qualitatively at higher resolutions. 



The computer calculations [^, agree quite well with experiment for different resolution parameters 
that implies the importance of the precise account for the energy-momentum conservation. 

Subjet multiplicities have also been studied for the separated quark and gluon jets. The analytical 
results ||149|| are seen (Fig. 12 from ||150|| ) to represent the data fairly well for large values of the subjet 
resolution scale Uq. 



3.10 Jet universality 

According to QCD, jets produced in processes initiated by different colliding particles must be universal 
and depend only on their own parent (gluon, light or heavy quark). This prediction has been confirmed 
by many experiments. For example, this is clearly seen if the multiplicities in the fragmentation region 
in the Breit system of ep-collisions are compared with the corresponding multiplicities in e"'"e~-processes 
at high energies. It has been found ||151|| that they coincide as is required by the universality condi- 



demonstrated in Fig. 8. Therefore we will not discuss it at some length but just stress once again the 
importance of such a prediction and its non-trivial origin. At the same time, let us mention that the 
universality of the properties of the jets produced does not imply the universality of the mechanisms of 
their production in various processes. 

4 Conclusions and outlook 

A list of successful analytical QCD predictions can be made longer. As was demonstrated above, quan- 
tum chromo dynamics has already predicted spectacular qualitative features of rather soft processes. 
Quantitatively, analytical results show that the higher order (NLO) terms always tend to improve the 
agreement with experiment compared to the asymptotic (LO) predictions. The accuracy achieved is 
often better than 20% or even 10% that is surprising by itself considering the rather large values of the 
expansion parameter of the perturbative approach. Moreover, some characteristics are very sensitive to 
ever higher order terms and should be carefully studied. The astonishing success of the computer so- 
lutions and Monte Carlo schemes demonstrates the importance of using the correct borderline between 
the perturbative and non-perturbative regions, which is approximately accounted in the analytical per- 
turbative approach by the cut-off parameter Qq and by the limits of integration over parton splitting 
variables. One can expect that the purely perturbative description becomes dual to the sum over all 
possible hadronic excitations. Nevertheless, the correspondence between the parton stage of the cas- 
cade evolution and the hadronization is not always defined at the quantitative level and requires further 
studies and the development of the common point of view. In particular, this problem is approached 
in a different way not only in analytical results if they are compared with the Monte Carlo models' 
conclusions but also in different Monte Carlo schemes. 

A new era of multiparticle production studies opens with the advent of new accelerators RHIC, 
LHC, TESLA, NLC, JLC, CLIC. We come closer to the asymptotic region^ even though the approach 
to the asymptotic laws will be, probably, extremely slow because all the predicted energy dependences 
of the physically measurable quantities are slow as well. Nevertheless, some predictions differ for various 
analytical approaches and Monte Carlo schemes at these energies and will be confronted to experimental 
data. It will allow to distinguish between them. The qualitative QCD predictions indicate the tendencies 
towards the asymptotic region where the perturbative estimates become more precise. 

The mean multiplicities will increase drastically. Now, in Au-Au collisions at the center of mass 
energies of 130 GeV per nucleon at RHIC the mean charged multiplicity exceeds 4000. It implies that 
the event-by-event analysis of various patterns formed by particles in the available phase space becomes 
meaningful. Such an analysis would allow one to classify in more detail the multiparticle production 
processes than, say, it is done in hadronic interactions where simple separation of the diffractive and 
non-diffractive processes is only considered. The study of the topology of individual events must provide 
us with a much richer information on the dynamics of the process compared with the results of the 
measurement of characteristics averaged over the whole event sample. Its results can be compared 
to the exclusive probabilistic Monte Carlo schemes. Searches for the super symmetric partners of well 
known particles, Higgses, new states of matter, new collective and interference effects, physics of jets 
and mini-jets will, surely, be among the most important directions of further investigations. 

The event-by-event approach would allow one to analyze the small colour-suppressed effects, which 
show the difference between the perturbative QCD and Monte Carlo calculations for the topology of 

"'^'^Let allow me to make here some short historical and lyrical digression. About 45 years ago, the tuitor of my diploma 
work I.Ya. Pomeranchuk told me that with the advent of the 10 GeV Dubna accelerator we enter the asymptotic region 
because 10 S> 1 (By 1 he meant the nucleon mass.). Now, we understand that this hope was too naive. Probably, our 
today's feelings are not more satisfactory even though they are supported by the modern QCD predictions. However, 
possible effects of the high parton density and of the non- linear interactions, which nowadays escape the detailed treatment, 
could change the situation at high energies. 



the individual events, properties of minijets or clusters (with the attraction-repulsion transition), other 
collective effects like the elliptic flow (and even the higher Fourier expansion terms of the azimuthal 
distribution in an individual event), the possible azimuthal asymmetry of the opening angle for indi- 
vidual jets, the ring-like events (the probable signature of the confinement and/or of the "Cherenkov 
gluons"), further analysis of the difference in accompanying gluon emission by heavy and light quarks 
(in particular, the angular distribution difference) etc. 

Of the principal importance is the study of the multiparticle correlations. The use of traditional 
formulas for the correlation functions is hampered by the large number of the independent variables. 
One can overcome this difficulty in spite of the seemingly complicated structure of these functions. 
The event-by-event analysis of experimental exclusive data can become available, quantified locally 
and provide statistically significant results if one uses wavelets (sometimes called the "mathematical 
microscope") for the pattern recognition in individual events ||152|| . With the help of wavelets one can 



separate correlations at the different resolution levels (from the short-range to long-range correlations) 
locally and in a compact form. In principle, the distributions of the wavelet coefficients can replace the 
complicated expressions for the multiparticle correlators. 

At the same time, the intriguing data obtained in the cosmic ray studies (for the recent review 
see Ref. | |153| ]) tell us that one can await for new phenomena in the fragmentation region of hadronic 
(and nuclear) collisions which, unfortunately, has not been devoted sufficient attention in accelerator 
experiments because most detectors are suited for the central range of rapidities. In particular, it would 
be important to learn the ratio between the baryon, meson and photon contents of the multiparticle 
production events. 

To confront QCD predictions with new experimental findings at ever higher energies will be crucial 
for the search for possible new states of matter. This is also important for the separation of any new 
physical signals from the conventional QCD background. 
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Figure captions 



Fig. 1. The energy dependence of average multiplicity of charged particles in e^e~-annihilation. 
The results of different fits according to formulas of perturbative QCD and of the Monte Carlo models 
are shown ( the solid and dotted lines are the fits of formula ( |23| ) with one and two adjusted parameters, 
the dashed line is given by the HERWIG Monte Carlo model; the vertically shaded area indicates the 
gluon jet data multiplied by the theoretical value of the ratio r (|2T|)). 

Fig. 2. The experimentally measured ratio r of multiplicities in gluon and quark jets as a function of 
energy in comparison with the predictions of analytical QCD and of the Monte Carlo model HERWIG 
(different QCD approximations, described in this paper, as well as 3NL0(e) with integration limits 
and l-e~^ in Eqns (^, (^ are indicated at the corresponding lines). 

Fig. 3. The ratio of the slopes of the energy dependences of mean multiplicities in gluon and quark 
jets according to experimental data and some theoretical calculations. 

Fig. 4. The ratio of the curvatures of the energy dependences of mean multiplicities in gluon and 
quark jets according to experimental data and some theoretical calculations. 

Fig. 5. The measured ratio Hq of the cumulant and factorial moments oscillates as a function of 
the rank q according to experimental data on multiplicity distributions of charged particles in e^e~- 
annihilation at the energy (the inset in the upper right corner shows the data for the moments of 
the ranks 2, 3 and 4). 

Fig. 6. The inclusive rapidity distribution of secondary particles in e+e~-annihilation has the shape 
of the hump-backed plateau (it is shown in the center of mass system). 



Fig. 7. The peak in the variable ^ is fitted at different energies by the distorted Gaussian with 
the moments predicted by the NLO-approximation of perturbative QCD (the fitted lines are drawn at 
A = 210 ) according to Q). 



Fig. 8. The peak position of the experimental inclusive distribution as a function of the mass of 
two jets is compared with NLO predictions (the central line is fitted by the CDF Collaboration data 
only). The lowest order (LO) predictions are shown by the lower straight line. The expected behaviour 
without the colour coherence effect gives rise to the upper line. 



Fig. 9. The energy dependence of the difference between the average multiplicities of charged parti- 
cles in events initiated by b and {u, d, s)-quarks. Experimental dots are compared with QCD predictions 
(the horizontal stripe) and with the results of the naive rescaling model (the decreasing stripe). 

Fig. 10. The charged hadronic flows in 3-jet events (the histogram) in comparison with analytical 
QCD predictions (the solid line) as functions of the azimuthal angle. 



Fig. 11. The normalized factorial moments of various ranks as functions of their scale (the size of the 
phase space window diminishes to the right on the abscissa axis). Different analytical approximations 
are compared with experimental data. LO-approximation (DLLA): (a) - ||128|| , (b) - [|130|| , (c) - ||129|| ; 
NLO-approximation (MLLA) - ||128|| . The qualitative but not quantitaive agreement is seen. 



Fig. 12. The subjet multiplicities in the separated gluon (a) and quark (b) jets as functions of the 
resolution parameter i/q are compared with analytical QCD results and with predictions of the Monte 
Carlo model JETSET for hadrons (HL) and partons (PL). 



